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Abstract. Recent results of particle production in the energy regime of the CERN-SPS 

are reviewed. In order to collect information on the properties of the QCD phase diagram 
^> ' systematic studies of the system size and the energy dependence of particle production in 

D t heavy ion collisions have been performed. Net-baryon distributions and results on strangeness 

production are discussed. The system size dependence of many observables can be understood 
O ' in the core-corona approach, which has implications on the possibility to use system size as a 

control parameter to study different areas of the phase diagram. Recent attempts to search 
S / for a critical point, such as multiplicity fluctuations and the transverse mass dependence of 

anti-baryon/baryon ratios are reviewed. 

>' 

^H . 

1. Introduction 

t^! ■ 

The QCD phase diagram, as it is illustrated in Fig.[JJ contains a variety of theoretically predicted 
features. The most important one is the phase boundary separating a hadron gas and a quark 
gluon plasma (red dashed line in Fig. [1]). Its position at vanishing baryonic chemical potential 
(fi B = 0) is defined by the critical temperature To. Even though recent attempts to calculate 
the exact value of Tq with lattice QCD still result in different numerical values [H [21 [3], there 
is a general consensus that for /i B = this transition is of the type of crossover. However, this 
crossover line might turn for a certain fi B into a first order phase transition, which would result 
in a critical point at that position [4J. However, predictions on the location of the critical point 
are still very uncertain [5]. It might very well be that it does not exist at all [6]. Recently, 
it has been conjectured that the phase diagram might contain yet another phase, the so-called 
quarkyonic phase [7], which would be separated from normal hadronic matter by a different 
phase boundary line (blue dashed line in Fig. [1]) . 

Experimentally, there are several control parameters available that might allow to study 
different regions of the phase diagram with heavy ion collisions. One is the variation of the 
center-of-mass energy which will force the reaction systems to follow different trajectories in the 
T-jUb plane, reflected by a change of the chemical freeze-out parameters. The second control 
parameter is the variation of the system size, which can either be achieved by performing central 
collisions of nuclei of different size, or by studying centrality selected minimum bias collisions. 
However, whether this can be considered as a good control parameter for probing different areas 
of the phase diagram shall be discussed further below. 
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Figure 1. The phase diagram of strongly interacting matter. The picture is taken from [8]. 



In the following we review some results from systematic studies of particle production in 
heavy ion collisions performed at the CERN-SPS, and will discuss what can be learned from 
existing data on the experimental accessibility of different areas in the phase diagram. Finally, 
we also present some results from existing attempts to search for a critical point. 

2. Variation of Energy 



Recent results on rapidity spectra of (anti-)protons in central Pb+Pb reactions, measured at 
energies between 20 A - 158 A GeV at the SPS [SI [TO], together with central Au+Au data from the 
AGS [11J and RHIC [TO] allow to study the energy evolution of stopping. Based on the measured 
rapidity spectra for p, p, A, A, H — , and H + , all corrected for feed down from weak decays, the 
net-baryon distributions B — B have been constructed [13]. The contribution of unmeasured 
baryons (n, S 1 * 1 , H°) was estimated using the results of a statistical hadron gas model |14j . In 
the SPS energy region a clear evolution of the shape can be observed (see left panel of Fig. [2]), 
changing from a single maximum at midrapidity towards a structure with two maxima and a 
dip at y = 0. Concurring with this evolution, also a strong change of the anti-baryon/baryon 
ratios takes place in the energy region of the SPS [101 [TO] and thus of the baryonic chemical 
potential at freeze-out /x B . 

The right panel of Fig. [2] shows a comparison of the energy dependence of mid-rapidity A, A, 
E~ , and H + production to several models. While the transport models UrQMD1.3 [16J and HSD 
|17] provide a reasonable description of the A/tt and A/ir ratios, they are clearly below the data 
points in case of the H~ and H + . This might indicate that an additional partonic contribution is 
necessary to reach the production rates observed for multi-strange particles. Statistical models 
on the other hand, generally provide a better match to the data. These models are based on 
the assumption that the particle yields correspond to their chemical equilibrium value and can 
thus be described by the parameters temperature T, baryonic chemical potential /j, B , volume 
V, and, in some implementations, by an additional strangeness under-saturation factor 7 S . The 
curves shown in the right panel of Fig. [21 labeled SHM(B), are taken from [18J and are based on 



parametrizations of the 



dependence of T and fx B , assuming 7 S = 1. The results of various 
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Figure 2. Left: The rapidity distributions of net-baryons at SPS energies |13|, I5| HP], together 
with results from the AGS [TT] and from RHIC [T2] for central Pb+Pb(Au+Au) collisions. 
Right: The energy dependence of the rapidity densities dn/dy around mid-rapidity of A, A, H - , 
and H + , divided by the total pion rapidity densities (n = 1.5 (ir~ + vr + )) for central Pb+Pb and 
Au+ Au collisions [IS1 |2D1 I2H 1221 HS1 [221 |2S I2S1 I2S1 I2Z1 1251 |2S1 ISS3I3I]- Also shown are results 
for the transport models UrQMD1.3 (dashed line) [TBJ and HSD (dotted line) [17], as well as a 
statistical hadron gas model (solid line) |18j . 



analysis of the measured particle abundances with statistical model approaches show that T and 
/x B parameter for the SPS energy range vary over a wide region, following an universal freeze-out 
curve [32], and might be in the vicinity of the possible critical point location [181 E3 EH ES]. 
Thus, the variation of J 
T-/i B plane. 



provides a well defined way of selecting different regions in the 



3. Variation of System Size 



3.1. Centrality Dependence of Net Protons 

Recent data on proton and antiproton production in minimum bias Pb+Pb reactions allow 
to study the system size dependence of stopping. Figure [3] shows the net-proton rapidity 
distributions for five different centrality classes, selected from minimum bias Pb+Pb interactions 
at 40 A (left panel) and 158 A GeV (right panel) [36J. A remarkable feature of this data is 
that there is no drastic change with centrality of the shapes of the distributions inside the 
measured region at 158^4 GeV. However, at 40^4 GeV the shapes of the distributions change 
with centrality, going from a two-maxima structure with a narrow dip at midrapidity for central 
collisions to a shallow valley like distribution for peripheral collisions, similar as observed at the 
AGS [37]. 

The p/p ratio, being directly sensitive to the baryonic chemical potential fi B , depends much 
less on system size than it does depend on energy. While for the mid-rapidity ratios at 158^4 GeV 
only a change in the order of a factor of two is observed when comparing central Pb+Pb reactions 
(p/p = 0.56, - 5% most central) to very peripheral ones (p/p = 0.11 for 43 - 100%), this ratio 
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Figure 3. Net-proton rapidity distributions for Pb+Pb collisions at 40A (left) and 158^4 GeV 
(right) for different centralities as measured by the NA49 experiment [36] • The full symbols 
represent the measurements. The open symbols are reflected at midrapidity. 



varies by a factor of 43 between central Pb+Pb collisions at 20^1 GeV (p/p = 0.0013, - 7%) 
and 158A GeV [TO]. Thus, only a small variation of // B can be achieved when studying reactions 
of different system sizes. 



3.2. Core Corona Model 

It has been found that the core-corona approach is able to describe the system size dependence 
of particle production reasonable well [381 EH HQl ED H2]- In this model a nucleus-nucleus 
collision is decomposed into a central core, which corresponds to the large fireball produced 
in central A+A collisions, and a peripheral corona, that is equivalent to independent nucleon- 
nucleon reactions. To quantify the relative contribution of the two components, the fraction of 
nucleons that scatter more than once /(A part ) can be used. f(N paIt ) can simply be calculated 
with a Glauber model [15]. This quantity allows for a natural interpolation between the yields 
Y measured in elementary p+p (= Icorona) and in central nucleus-nucleus collisions (= lcore) : 



Y(N paxt ) = ^Part [/(A part ) >core + (1 - /(A part )) ^corona] 



1 



The left panel of Fig. [4] shows for example the result of this approach compared to the system 
size dependence of the strangeness enhancement factors E at 158A GeV |44[ 145] . Here, E is 
defined relative to p+p reactions as baseline measurement: 
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The observed rapid rise for very small systems and the subsequent saturation can naturally be 
described within this model. 

It is interesting to observe that this approach not only works for yields, but also for dynamical 
quantities such as the average transverse momenta (m t ) — wiq (see right panel of Fig. [4|). This 
suggests that the core-corona approach in general provides a reasonable way for understanding 
the evolution from elementary p+p to central Pb+Pb collisions. 

A more differential look on the core-corona model is shown in the left panel of Fig. [5] for 
symmetric reaction systems. Due to their different surface to volume ratio, the A par t dependence 
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Figure 4. Left: The midrapidity yields per wounded nucleoli relative to p+p yields for central 
C+C, Si+Si and minimum bias Pb+Pb reactions at 158^4 GeV [441 I45j . Right: The {nit) — too 
values at mid-rapidity for Pb+Pb collisions at 40^4 and 158A GeV, as well as for near-central 
C+C and Si+Si reactions at 158^4 GeV [3]. The (anti-)proton data are taken from [46j. Also 
shown are the results from a fit for A and A with the core-corona approach (solid lines). 



of /(A r part ) is much steeper for smaller reaction systems than for larger oneaj Also, the maximum 
value of /(Ap art ) depends on the size of the nuclei. E.g., while for very central Pb+Pb collisions 
/max(Ap art ) ~ 0.9 can be reached, the maximum value for very central C+C reactions is 
/max(Ap art ) « 0.65. This A dependence of f max is illustrated by the dashed red line in the left 
panel of Fig. [SJ This has the consequence that the different relative core corona contributions 
have to be taken into account when comparing even very central nucleus-nucleus collisions of 
different size. Especially, this might explain any apparent system size dependence of the chemical 
freeze-out parameters T and [X B , as determined by statistical model fits to central A+A collisions 
of different size (e.g. |14|). Therefore, we conclude that, assuming the core-corona model is valid, 
the system size might not provide a good control parameter to probe different regions of the 
QCD phase diagram with produced fireballs of different temperature. It is more likely that one 
only observes a change in the relative admixture of central fireball (core) and peripheral p+p like 
corona, whose freeze-out parameters are different, but independent of the size of the involved 
nuclei. 

Another interesting aspect of the core-corona model is shown in the right panel of Fig. [5j 
Here, the A part dependence of /(A part ) is shown for asymmetric collisions. In this case quite 
distinct centrality dependences can be observed. While for symmetric collisions this dependence 



1 Please note that the curves shown in Fig. f4]are based on a function /(iV par t) that was calculated for Pb+Pb 
interactions [JT]. Therefore their direct comparison to the smaller systems C+C and Si+Si is not possible. 
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Figure 5. The fraction f(N part ) of participating nucleons that scatter more than once as a 
function of the number of participants N paIi . f(N pa _ Tt ) was calculated with a Glauber model 
|43t I47] , The left panel shows results for symmetric systems, the right panel for asymmetric 
ones. The dashed lines connect the values for the most central collisions. 

is following a continous rise (left panel of Fig. [5]) , for asymmetric collision systems with a small 
projectile nucleus (e.g. O+Pb or Si+Pb) a rapid rise followed by a maximum and a subsequent 
decrease of /(iV part ) is seen. Therefore, a similar centrality dependence of particle yields can be 
expected for these type of collisions and its measurement would constitute a test for the validity 
of the core-corona model. 

4. Search for the Critical Point 



4-1. Fluctuations 

One of the observables proposed to be sensitive to the presence of a critical point in the QCD 
phase diagram are fluctuations of either mean transverse momentum or multiplicity [38]. If 
the chemical freeze-out is happening in the vicinity of the critical point, an enhancement of 
fluctuations might be visible. The NA49 experiment has therefore performed a systematic study 
of multiplicity fluctuations in very central (1% most central) Pb+Pb collisions as a function of 
beam energy |49j . Figure [6] shows the resulting scaled variance uj = Var(n)/(n) as a function of 
baryonic chemical potential [i B , that has been derived from statistical model fits to the particle 
abundances measured at the different beam energies [14J. 

While for all charged particles taken together u is close to unity, it is slightly below one for 
negatively and positively charged particles analyzed separately, i.e. the distributions are in these 
cases a bit narrower than the corresponding Poissonian. However, in none of the studied charge 
combinations a significant energy dependence has been so far observed. 

Figure [6] also shows some theoretical expectations for u>, which result from the combination 
of several assumptions [50]. The position of the critical point in terms of T and /i B is taken 
from a lattice QCD calculation [35]. The magnitude of the fluctuations at the critical point is 
based on |48t 151] , assuming two different correlation lengths £ = 3 fm (solid line) and £ = 6 f m 
(dashed line), while the widths of the enhancements around the critical point has been chosen 
according to [52J as cr(/j, B ) ~ 30 MeV. Both predictions are at variance with the data, so that 
up to now there is no evidence for a critical point in fluctuation measurements. 
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Figure 6. Energy dependence of multiplicity fluctuations, given by the scaled variance u>, 
for the 1% most central Pb+Pb collisions in the forward rapidity region (1.1 < y w < ybeam) 
as measured by the NA49 experiment |49t I5UJ. The different beam energies are represented by 
the corresponding fi B values, determined by a statistical model fit [14J. The lines correspond to 
predictions for for critical point (see text). 



4-2. Transverse Mass Spectra of Baryons and Anti-baryons 

Another possible way of detecting a critical point has been suggested in [53] . Here it was shown 
that the proximity of a critical point, or rather a critical area, might deform the isentropic 
trajectories of an expanding fireball in the T-fi B plane, due to a focussing effect towards the 
critical point. As a consequence one should see a decrease of the anti-baryon/baryon ratios with 
increasing transverse mass, while in the absence of a critical point this ratio should rather be 
flat as a function of m% . 

The left panel of Fig. [7J shows the p/p ratio versus mt — mo for five different beam energies 
|1U| I5l)| . In order to quantify any change in transverse mass dependence, a linear fit has been 
applied to the data. The A/ A and H + /H~ ratios have been analyzed in the same fashion. The 
resulting energy dependence of the slope parameters a is summarized in the right panel of Fig. [71 
No significant energy dependence can be seen and thus this observable currently does not provide 
any evidence for a critical point. 



5. Conclusions 



In order to probe different areas in the QCD phase diagram, the variation of the center-of-mass 
energy can be considered as the best defined control parameter. Especially in the energy regime 
of the SPS a wide region in the T-/i B plane can be explored experimentally, as is evident from 
the strong dependence of the chemical freeze-out parameters on the energy here. The situation 
is much less well defined when using the system size as control parameter. Since the system size 
dependences of many observables (e.g. strangeness production and (m t ) — mo) seem to be well 
described within the core-corona approach, it is unlikely that the freeze-out temperatures of the 
fireballs vary with their size. Rather one observes a change of the relative admixture of the core 
and corona contributions. However, the trajectories of the fireballs in the core might be always 
the same, irrespective of their size. Studying centrality dependences in asymmetric reaction 
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Figure 7. Left: The p/p ratio versus the transverse mass mt — mo |10j I50j . The solid lines 
show a linear fit to the ratios. Right: The energy dependence of the slope parameters a from 
the fits for different B/B ratios [TO1 [13 [50] . The insert shows the same data with an expanded 
vertical scale. 



systems could provide an additional experimental test of the core-corona model. In these cases 
this model predicts for smaller projectile nuclei a maximum in the relative core contribution for 
semi-central reactions, in contrast to symmetric reaction systems that reach the maximum for 
very central collisions. 

Several attempts have been made to search within existing data for evidences of a critical 
point, such e.g. as fluctuations and the transverse mass dependence of anti-baryon/baryon 
ratios. So far non of these searches yielded a positive result. The success of future searches 
will on one side mainly depend on the development of more refined and differential (i.e. scale 
and pt dependent) observables, that allow to distinguish effects due to a critical point from 
more trivial sources. On the other side, since the magnitude of these effects might be small, 
experimental progress on reducing systematic and statistical errors will be necessary. 
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